TAIL ESTIMATES FOR THE BROWNIAN EXCURSION 
AREA AND OTHER BROWNIAN AREAS 



SVANTE JANSON AND GUY LOUCHARD 

Abstract. Several Brownian areas are considered in this paper: the 
Brownian excursion area, the Brownian bridge area, the Brownian mo- 
tion area, the Brownian meander area, the Brownian double meander 
area, the positive part of Brownian bridge area, the positive part of 
Brownian motion area. We are interested in the asymptotics of the right 
tail of their density function. Inverting a double Laplace transform, we 
can derive, in a mechanical way, all terms of an asymptotic expansion. 
We illustrate our technique with the computation of the first four terms. 
We also obtain asymptotics for the right tail of the distribution function 
and for the moments. Our main tool is the two-dimensional saddle point 
method. 



1. Introduction 

Let i?ex(i)) t £ [0)1]) ^ (normalized) Brownian excursion, and let 
■Sex := /o Bex{t)dt be its area (integral). This random variable has been 



studied by several authors, including Louchard [13; UM, Takacs [19|], and 
Flajolet and Louchard [3]; see also the survey by Janson [§] with many 
further references. 

It is known that Bex has a density function /ex, which was given explicitly 
by Takacs (l9l ] as a convergent series involving a confluent hypergeometric 
function. (The existence and continuity of /ex follows also from Theorem [XT] 
below.) The series expansion of /ex readily yields asymptotics of the left tail 
of the distribution, i.e., of fex{x)aiid P(Sex < x) as x ^ O^see Louchard 
14 1 and Flajolet and Louchard [7| (with typos corrected in 9]). 



The main purpose of this paper is to give corresponding asymptotics for 
the right tail of the distribution of i3ex, i-e., for the density function fexix) 
and the tail probabilities P(^cx > x) for large x. We have the following 
result. 

Theorem 1.1. For the Brownian excursion area, as x —> cxd, 

/ex(x) ~ ^x^e-e-^ (1.1) 
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¥{Be.>x)^^xe-'^\ (1.2) 



and 



More precisely, there exist asymptotic expansions in powers of x ^, to arbi- 
trary order N , as x ^ oo, 

72\/6 2 -fi.T2 A 1 -2 5 _4 25 _Q . _2Af 



^ ^ V 9 1296 46656 ^ 

6\/6 _fi^2 1 „2 1 -4 7 -6 



P(B„ >.) = -I^.e-'.' ^1 + + O 

Unlike the left tail, it seems difficult to obtain such results from Takacs's 
formula for /ex, and we will instead use a method by Tolmatz [2^; [24; 25] that 



he used to obtain corresponding asymptotics for three other Brownian areas, 
viz., the integral ^Bbr := /q^ |^br(i)| dt of the absolute value of a Brownian 
bridge B\^j.{t), the integral B]^^ := \B{t)\dt of the absolute value of a 

Brownian motion B{t) over [0, 1], and the integral ;Bbr+ := Jq ^br(i)+ dt of 
the positive part of a Brownian bridge. 

The (much weaker) fact that — lnP(;Sex > x) ~ — 6x^, i.e., that 

P(5ex > x)= exp(-63;2 + o{x^)) , (1.3) 

was shown by Csorgo, Shi and Yor Q| as a consequence of the asymptotics 
of the moments E B^^ found by Takacs [l^ , see Section [9l It seems difficult 
to obtain more precise tail asymptotics from the moment asymptotics. It is, 
however, easy to go in the opposite direction and obtain moment asymptotics 



from the tail asymptotics above, as was done by Tolmatz [2J; |25| for Bi,^, 
i3br and ^br+; see again Section [9l In particular, this made it possible to 
guess the asymptotic formula (jl.ip before we could prove it, by matching 
the resulting moment asymptotics with the known result by Takacs [19]. 

An alternative way to obtain (II. 3p is by large deviation theory, which 
easily gives (II. 3p and explains the constant 6 as the result of an optimization 
problem, see Fill and Janson [£]. This method applies to the other Brownian 
areas in this paper too, and explains the different constants in the exponents 
below, but, again, it seems difficult to obtain more precise results by this 
approach. 

Besides the Brownian excursion area and the three areas studied by Tol- 
matz, his method applies also to three further Brownian areas: the integrals 

Bme ■= Jo \Bmeit)\dt, Bdm ■= /q^ ] -Bdm (*) ] dt and Bhm+ ■= /q^ -B(t)+ dt of 
a Brownian meander B^e{t), a Brownian double meander Bd^{t), and the 
positive part of a Brownian motion over [0,1]. We define here the Brow- 
nian double meander by -Bdm(i) := B{t) — mino<u<i i?(u); this is a non- 
negative continuous stochastic process on [0, 1] that a.s. is at a unique 
point T £ [0, 1], and it can be regarded as two Brownian meanders on the 
intervals [0, r] and [r, 1] joined back to back (with the first one reversed), see 
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Majumdar and Comtet [3] and Janson [§]; the other processes considered 
here are well-known, see for example Revuz and Yor 18l |. 

We find it illuminating to study all seven Brownian areas together, and 
we will therefore formulate our proof in a general form that applies to all 
seven areas. As a result we obtain also the following results, where we for 
completeness repeat Tolmatz's results. (We also extend them, since Tolmatz 
13; 24; IH] gives only the leading terms, but he points out that higher order 



terms can be obtained in the same way.) We give the four first terms in 
the asymptotic expansions; they can (in principle, at least) be continued to 
any desired number of terms by the methods in Section [6l only even powers 
x~'^^ appear in the expansions. 

Theorem 1.2 (Tolmatz [23]). For the Brownian bridge area, as x — > oo, 
2\/6 _6^2 A , 1 -2 , 1 ^-4 , 25 e , 



■'^'^ 0r V 18 432 46656 ^ 

P(^br > x) = ^ f 1 - 1 + ^ - + O (x-')] 

^ ' ^ V 36 108 46656 ^ ' ) 

Theorem 1.3 (Tolmatz [24]). For the Brownian motion area, as x ^ oo. 

Theorem 1.4. For the Brownian meander area, as x ^ oo, 

= 3^/3.e--'/. (l - 1 X- - + ^ + O (.-)) . 

F(fi„. > X) = V3e--'/. (i _ ^ + _! _ H + o (.-)) . 
Theorem 1.5. For the Brownian double meander area, as x ^ oo. 



F,Ba„ > .) ^ ^ .-e-V= (, _ 1 + ^ _ |1 + o (.-)) , 
Theorem 1.6 (Tolmatz [is']). For ^/le positive part of Brownian bridge area. 



as X ^ oo, 



^ -6x2 / 1^-2 7 4 17 e . 8, 



2^/6¥ V 18 5184 186624 
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Theorem 1.7. For the positive part of Brownian motion area, as x ^ oo, 

^ ^ e-V. (l + i - A + ^ + O (.-)) 

It is not surprising that the tails are roughly Gaussian, with a decay like 

2 

g-cx £qj. gQ]-Qg constants c. Note that the constant in the exponent is 6 for 
the Brownian bridge and excursion, which are tied to at both endpoints, 
and only 3/2 for the Brownian motion, meander and double meander, which 
are tied to only at one point. It is intuitively clear that the probability of 
a very large value is smaller in the former cases. There are also differences 
in factors of x between B^r and Bex, and between Bhm and Bme, where the 
process conditioned to be positive has somewhat higher probabilities of large 
areas. These differences are in the expected direction, but we see no intuitive 
reason for the powers in the theorems. We have even less explanations for 
the constant factors in the estimates. 

Remark 1.8. If we define ^br- := -Bbr(t)- dt, we have ;Bbr = 'Bbr++'Bbr-; 
further, B^r- = Bhv+ by symmetry. Hence, for any x, 

F{Bhr > x)> P(Sbr+ > X or Bhr- > x) 

= F{Bbr+ >X)+ P(Sbr- > X) - P(Sbr+ > X and Sbr- > x) 

> 2F{Bbr+ >x)- 2F{Bbr > 2x). 

By Theorems 11.21 and 11.61 the ratio between the two sides is 1 + ^x^'^ + 
0{x~^); hence, these inequalities are tight for large x. This shows, in a very 
precise way, the intuitive fact that the most probable way to obtain a large 
value of ;Sbr is with one of ^br+ and ^br- large and the other close to 0. 

The same is true for ^Bbm and ;Bbm± by Theorems 11.21 and 11.61 It is 
interesting to note that for both ^Bbr and ^bmi the ratio F{B > x)/2F{B+) = 
1 + ^x^'^ + 0{x~^), with the first two terms equal for the two cases (the 
third terms differ). 

Tolmatz's method is based on inverting a double Laplace transform; this 
double Laplace transform has simple explicit forms (involving the Airy func- 
tion) for all seven Brownian areas, see the survey [9] and the references given 
there. The inversion is far from trivial; a straightforward inversion leads to 
a double integral that is not even absolutely convergent, and not easy to 
estimate. Tolmatz found a clever change of contour that together with prop- 
erties of the Airy function leading to near cancellations makes it possible to 
rewrite the integral as a double integral of a rapidly decreasing function, for 
which the saddle point method can be applied. (Kearney, Majumdar and 



Martin [I4] have recently used a similar change of contour together with 



similar near cancellations to invert a (single) Laplace transform for another 
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type of Brownian area.) We follow Tolmatz's approach, and state his in- 
version using a change of contour in a rather general form in Section [3l the 
proof is given in Section [8l This inversion formula is then applied to the 
seven Brownian areas in Sections HHSl Moment asymptotics are derived in 
Section [9l 

A completely different proof for the asymptotics of F(Bhr > x) and 
P(^bm > x) in Theorems 11.21 and 11.31 has been given by Fatalov [5] using 
Laplace's method in Banach spaces. This method seems to be an interest- 
ing and flexible alternative way to obtain at least first order asymptotics in 
many situations, and it would be interesting to extend it to cover all cases 
treated here. 

We use Ci, C2, . . . and ci, C2, . . . to denote various positive constants; ex- 
plicit values could be given but are unimportant. We also write, for example, 
Ci(M) to denote dependency on a parameter (but not on anything else). 

2. Asymptotics of density and distribution functions 

The relation between the asymptotics for density functions and distribu- 
tion functions in Theorems I1.1H1.7I can be obtained as follows. 

Suppose that X is a positive random variable with a density function / 
satisfying 

/(x) ~ax"e-^^', x^oo, (2.1) 

for some numbers a,b > 0, a £ M. It is easily seen, e.g. by integration by 
parts, that ()2.ip implies 

P(X > x) ~ ^x^-^e"^^', x^oo. (2.2) 

Obviously, there is no implication in the opposite direction; X may even 
satisfy (j2.2p without having a density at all. On the other hand, if it is 
known that (j2.ip holds with some unknown constants a, b, a, then the 
constants can be found from the asysmptotics of ¥{X > x) by (|2.2p . 

The argument extends to asymptotic expansions with higher order terms. 
If, as for the Brownian areas studied in this paper, there is an asymptotic 
expansion 

fix) = x"e-^^' (ao + a2X-2 + qax''^ + ■■■ + 0{x-^^)) , x ^ 00, (2.3) 

then repeated integrations by parts yield a corresponding expansion 

F{X > x) = x^-^e-^^" (a'o + a'^x'^ + a'^x'^ + ■■■ + 0{x-^^)) , x ^ 00, 

(2.4) 

where g'q = ao/(25), a'2 = ao{a — l)/(25)^ -|- 02/(26), ...; in general, the 
expansion (j2.3p is recovered by formal differentiation of (j2.4p . which gives a 
simple method to find the coefficients in (12. 4p . 
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3. A DOUBLE Laplace inversion 

We state the main step in (our version of) Tolmatz' method as the fohow- 
ing inversion formula, which is based on and generaUzes formulas in Tolmatz 



23; M,\m. 



oo 



ds = ^(x), x>0, (3.1) 



Fractional powers of complex numbers below are interpreted as the prin- 
cipal values, defined in C \ (— oo,0]. 

Theorem 3.1. Let X be a positive random variable and let il){s) := Ee"'^"''' 
be its Laplace transform. Suppose that < < 3/2 and that 

where ^ is an analytic function in the sector {z £ C : \ argzj < 57r/6} such 
that 

^{z) = odzl-"), z ^0 with] avgz\ <57r/6, (3.2) 

^'(z)=0(l), \z\ ^ oo with \argz\ <5Tr/6. (3.3) 

Let 

q,*{z) := e2^'^'/3^(e2^'/3^) - e-^'"''/^^{e-^'''/^z) . (3.4) 
Finally, assume that 

^r*(^) = 0(|^|-6)^ \z\ ^ OO with \aTgz\ < IT /6. (3.5) 

Then X is absolutely continuous with a continuous density function f given 
by, for X > and every ^ > 0, 

■ \ / expfe-2/3gecee'^-e^^(esec0)3/V3/2 

Je=-TT/2 Jr=Q ^ 

■ e(^-2^/3)'^(sec0)^/2-!.^i.-5/2^*(^gie/3)^^^^_ (3_g) 
Note that ^* is analytic in the sector |arg2:| < -I, with, by ()3.2p and 

^*{z) = o{\z\-''), z^O with largzl < -, (3.7) 

6 

■^*(z) = 0{l), Izl ^ oo with I argzl < -. (3.8) 

6 



However, we need, as assumed in (j3.5p . a more rapid decay as |z| ^ oo than 
this. 

Remark 3.2. In all our applications, ^' is, in fact, analytic in the slit plane 
C \ (— oo, 0], and (j3.2p and (j3.3p hold in any sector | arg z\ < it — 6; thus ^* 
is analytic in | argz| < tt/3, and (j3.7p and (|3.8p hold for | argzj < vr/S — S. 



Remark 3.3. To obtain Tolmatz' version of the formulas, for example [2S 
(30)] (correcting a typo there), take z/ = 1/2 and "if* as in (|4.5p below, and 
make the substitutions x = A, ^ = aX'^^^ and r = ap~'^^^ sec 6. 
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We prove Theorem 13.11 in Section [HI but show first how it applies to the 
Brownian areas. 

4. The function for Brownian areas 
For the Brownian bridge area i3br we have v = 1/2 and, see e.g. [9, (126)], 

= -2V6^i(^ (4 1) 

Ai'(2i/3^)' ^^-'^ 

which by the formula [1, 10.4.9] 

Ai(ze^2W3) = 1 g±W3 (Ai(z) ^ iBi(z)) (4.2) 
and its consequence 

Ai'(ze±2W3) = ie^"^/3(Ai'(^) ^ iBi'(z)) (4.3) 
together with the Wronskian III, 10.4.10] 

Ai(z)Bi'(z) - Ai'(z)Bi(z) = vr"^ (4.4) 

by a simple calculation leads to, as shown by Tolmatz [23|, Lemma 2.1], see 
(itni below, 

" Ai'(2i/3^)2 + Bi'(2V3^)2- (^•^) 

It seems simpler to instead consider y^^Bbr- Note that, by the simple change 
of variables s ^ and x i— > 2~^/^x in (13. 1|) . if (13. ip holds for some 

random variable X and a function it holds for \/2X and 2^'^/'^^(2^^/'^z). 
We use the notations ^'br and ^^^^ for the case X = \/^Bhr and obtain from 
(j4.ip the simpler 

= (4.6) 

and thus, by (jM]), (SiS]) and (|01) . 

^L(^) = =E^e^3W3 Ai(^)TiBi(^ 



^ ^ Ai'(z)TiBi'(z) 

^ (Ai(z) T iBi(z)) (Ai^(z) ± iBi^(z)) 
^ Ai'(z)2 + Bi'(z)2 



(4.7) 



Ai'(z)2 + Bi'(z)2' 



For the Brownian excursion area B^x we have = 1/2 and by Louchard 
il, see also @, (80)], 

Mz) - -2'/'- ( ^^i^] - 2V2 foi/3 Ai-(2V3.) \ ' _ 3/2^ (4 8) 
^ dz l^Ai(2V3^) j Ai(2V3^) J ^' ^^-^^ 
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Again, it seems simpler to instead consider ^/2B^^, for which we use the 
notation ^'ex and ^l^. We have, see Louchard [M] and [9., (81)], or by (f48]l 
and the general relation above, 

d /Ai'(z)\ fA.Hz)V , , 

and thus by ([32]), and ([OD 

± 

= y ± f e±-V32 f,T2W3Ai^(^)TiBi-(.)\ ^ _ 2^i2.i/3^ 

^ V V Ai(z) ^ iBi(z) ; 

= 2 y ^,T3.i/3 (Ai^(^) T iBi-(z))^(Ai(z) ± iBi(z))^ ^ ^ 
^ (Ai(z)2+Bi(z)2)2 

^2y- (Ai^(z)Ai(z) + Bi^(z)Bi(^)Ti^-^)' 

±^ (Ai(z)2+Bi(z)2)' 
_ 87r-M(Ai(z)Ai'(z) + Bi(z)Bi'(z)) 
(Ai(z)2 +Bi(z)2)2 



(4.10) 



The Brownian motion area i3bm is another case treated by Tolmatz 24 1. 
Note that in this case v = \. For \/2;Bbmi we have by Takacs [2?], see also 
Kac [To], Perman and Wellner [17j], and [S, Section 20 and Appendix C.l], 
V = \ and 

^bm(^) = -JJ^y (4-11) 
where we use the notation, see [ol, Appendix A], 

Al{z):= Ai{t)dt = -- Ai{t)dt. (4.12) 

Jz 3 Jo 



If we further define 



Bl{z) := I Bi{t)dt, (4.13) 
Jo 



we have by (j4.2p 

±2^i/3 

AI(ze±2W3) = - - / Ai(t)dt 
3 Jo 



l_e±2W3 rAi(te±2-/3)dt 
3 Jo 

g±3W3 r(Ai(t)^iBi(t))dt 

Jo 



3 

1 1 
3 ~ 2 

i-iAI(z)T^iBI(z). (4.14) 
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± 

^ V ^ +3W3 l-AI(^)TiBI(z) 
Ai'(z) ^ iBi'(z) 

1 - Al{z) T iBI(z)) (Ai'(z) ± iBi'(z)) 



Ai'(z)2 + Bi'(z)2 



^. Bi^(z)-AI(z)Bi^(z)-Ai^(z)BI(z) 

Ai'(z)2 + Bi'(z)2 ■ ^ ■ ^ 



For the Brownian meander, or more precisely \/2i3nie) by Takacs 22], see 
also [1, Section 22 and Appendix C.3], (|3.ip holds with u = 1/2 and 



^me(^) = (4.16) 

Ai(z) 



Consequently, using ()4.14p and ()4.2p . 

^^^l-AI(z)TiBI(z) 



^ Ai(z) ^ iBi(z) 

(1 - AI(z) ^ iBI(z)) (Ai(z) ± iBi(z)) 



^ Ai(z)2 + Bi(z)2 

^. Bi(z)-AI(z)Bi(z)-Ai(z)BI(z) 

Ai(z)2 + Bi(z)2 ■ ^ ' ' 



For the Brownian double meander, or more precisely ^/2Bdm, by Ma- 
jumdar and Comtet [1^, see also [i. Section 23], 1^ holds with = 1 
and 
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Consequently, using ()4.14p and (14. 2p . 

± 

^ V Ai(z)TiBi(z) ) 
^ ^ ^ ((l-A I(z) T iBI(z)) (Ai(z) ± iBi(z)))' 



± 



(Ai(z)2 + Bi(z)2)2 



. ((1 - AI(z))Ai(z) + BI(z)Bi(z)) ((1 - AI(z))Bi(z) - BI(z)Ai(z)) 
' (Ai(z)2 +Bi(z)2)2 ■ 

(4.19) 

The positive part of a Brownian bridge is another case treated by Tolmatz 



251]. For \/2;Sbr+, by Perman and Wellner 17 1, see also Tolmatz [25i | and 



Section 22 and Appendix C.2], (j3.ip holds with = 1/2 and 

Consequently, by ([321), (USD and (OI) . 

= V ^,,.^i/3 Ai(e^2.i/3,) 

e±7ri/32l/2Ai(e±2-i/32) - Ai'(e±2-i/3^) 



2E± 

± 
± 

2E± 



Ai(z) =F iBi(z) 



zV2(Ai(z) ^ iBi(z)) - eT3W3(Ai'(z) ^ iBi'(z)) 

Ai(z) =F iBi(z) 
(zi/2Ai(z) + Ai'(z)) ^ i(zi/2Bi(z) + Bi'(z)) 

(Ai(z) ^ iBi(z)) (zi/2Ai(z) + Ai'(z) ± i(zi/2Bi(z) + Bi'(z))) 



± (2V2Ai(z) + Ai'(z))^ + (zi/2Bi(z) + Bi'(z))^ 

4i7r-i 

ri/2Ai(z) + Ai'(z)) ^ + (zV2Bi(z) + Bi'(z)) ^ ■ (4-21) 



For the positive part of a Brownian motion, or more precisely 



bm+5 



by Perman and Wellner [17[, see also [S, Section 23 and Appendix C.l], (j3.1jl 
holds with V = 1 and 

_ z-V2Ai(z) + AI(z) 

- .V2Ai(.)-Ai'(.)- (^-^^^ 
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Note that this ^ is singular at 0, but still satisfies ([321) • By ([M]), <^M, 
(f4T4]l and g31), 



± 

_ V- ±2.i/3 e^-^/3^-i/2Ai(e±2M/3^) + Al(e±2-/3^) 
- Z^^^ e±W32i/2Ai(e±2W3^) _ Ai'(e±2W3^) 

_ ^ (Ai(z) ^ iBi(z)) + 1 - Al{z) T iBI(z) 

± 



zV2(Ai(z) ^ iBi(z)) - eT3W3(Ai'(z) If iBi'(z)) 

_ ^ z-i/2Ai(z) + 1 - AI(z) ^ i(z-V2Bi(z) + BI(z)) 
" ± ^ (zV2Ai(z) + Ai'(z)) ^ i(zi/2Bi(z) + Bi'(z)) 

. (z-V2Ai(z) + 1 - AI(z)) (zi/2Bi(z) + Bi'(z)) 
' (zi/2Ai(z) + Ai'(z))^ + (zi/2Bi(z) + Bi'(z))^ 

^. (z-i/2Bi(z) + BI(z)) (zi/2Ai(z) + Ai'(z)) 
' (zV2Ai(z) + Ai'{z)f + (zi/2Bi(z) + Bi'(z))^ 
^. (1 - AI(2)) (zi/2Bi(z) + Bi'(z)) - Bl{z) (zV2Ai(z) + Ai'(z)) + z-^/^-r-^ 



^i/2Ai(z) + Ai'(z))^ + (zi/2Bi(^) + Bi'(z))^ 



(4.23) 



Note that the functions ^'br, ^^ex, ^bm; ^me and ^'dm given above in (j4.6p . 
(j4.9p . (14. lip . (I4.16[) . (j4.18p are meromorphic, with poles only on the negative 
real axis, because the only zeros of Ai and Ai' are on the negative real axis 
[H, p. 450]. The functions ^'br+ and ^'bm+ in ([420]) and (14221) are analytic in 
the slit plane C \ {—oo, 0], since Tolmatz (25l ] showed that z^/^Ai(z) — Ai'(z) 
has no zeros in the slit plane; see Appendix |A] for an alternative proof. In 
particular, all seven functions are analytic in the slit plane. Furthermore, 
all except ^bm+ have finite limits as z — > 0, and in particular they are 0(1) 
as z ^ so 1^ holds. By (|02]) . we have ^'bm+(2) ~ 2"^/2Ai(0)/Ai'(0) 
and thus ^'bm+ = 0(|z|~^/^) as z — > 0; since in this case u = I, (|3.2p holds 
for ^'bm+ too- 
Next we consider asymtotics as \z\ oo. The Airy functions have 
well-known asymptotics, see [H, 10.4.59, 10.4.61, 10.4.63, 10.4.66, 10.4.82, 
10.4.84]. The leading terms are, as \z\ oo and uniformly in the indicated 
sectors for any 5 > 0, 

-1/2 

Ai(z) - l_2-i/4e-2^«/V3^ I arg(2)| <TT-6, (4.24) 

— 1/2 

Ai'(z) ~ _!L_^i/4g-2.^/V3^ I arg(z)| < vr - <5, (4.25) 

-1/2 

AI(z) ~ z-3/4e-2^3/V3^ I arg(z)| < vr - 5, (4.26) 
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Bi'iz 



Bl(z) ~ ^-1/2^-3/4^2.3/2/3 



arg(z)| < 7r/3 - 5, (4.27) 
arg(z)| < 7r/3 - 5, (4.28) 
arg(z)| < 7r/3 - 5. (4.29) 



It follows by using (|CTD . (ITOD and (IM]) in (gSD, (031), (gTH), 
(|4.20p . (j4.22p that in all seven cases (j3.3p holds; more precisely, ^'(2:) ~ z~'^ 
as l^l — > 00 with I arg2:| < n — 6. (For real 2; > 0, this is always true, as 
follows from (|3.ip by the change of variables s = t/x and monotone (or 
dominated) convergence.) 

Turning to ^* , we observe first that, by (j3.4p . in all seven cases, ^*{z) 
is analytic in | argz| < 1/3. Next, ()4.24p - (14.29p show that, as \z\ — > 00 in 
a sector |arg(z)| < 7r/3 — 5, Ai, Ai',AI decrease superexponentially while 
Bi, Bi',BI increase superexponentially. Hence, we can ignore all terms in- 
volving Ai. More precisely, gZD, dTO]) - (ITO]) . (HTTD . dMU, (03]) 
together with (j4.24p - (j4.29p yield the asymptotics, as \z\ — > 00 with (for 
example) | argzj < vr/G, 



^l,{z 

*bm+(^ 



27r-ii 
BiW 



i + ofe-«-^^V3)y 



(l + 0(e-«^^^V3)), 



Bi'(z 

21 
BiR 

4i 



l + 0(e-2^^^^/3)). 



Bi(z)2 

4i7r-i 



(zi/2Bi(z) + Bi'(z)) 
2i 



1 + o (e-^^"'l^ 



ji/2Bi(z) + Bi'(z) 



1 + e 



-2.3/2/3 



(4.30) 
(4.31) 
(4.32) 
(4.33) 
(4.34) 
(4.35) 
(4.36) 



In all seven cases, ^* decreases superexponentially in the sector; in par- 
ticular, (|3.5p holds. It is remarkable that in all seven cases, ^{z) decreases 
slowly, as ^-^/^ or z^^ , but the linear combination ^'*(z) decreases extremely 
rapidly in a sector around the positive real axis; there are thus almost com- 
plete cancellations between the values of ^{z) at, say, arg z = ib27ri/3. These 
cancellations are an important part of the success of Tolmatz's method. 

We have verified all the conditions of Theorem 13.11 Hence, the theorem 
shows that the variables have continuous density functions given by (|3.6p . 
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5. The saddle point method 



We proceed to show how the tail asymptotics for the Brownian areas 
follow from Theorem 13.11 and the formulae in Section |4] by straightforward 
applications of the saddle point method. For simplicity, we give first a 
derivation of the leading terms. In the next section we show how the cal- 
culations can be refined to obtain the asymptotic expansions in Theorems 
[IIHLTl 

We use H G {br, ex, bm, me, dm, br+, bm+} as a variable indicating the 
different Brownian areas we consider. We begin by writing (j4.30p - (j4.36p . 
using (|i:27|l and (fOHl) . as 



(5.1) 



where 7br = 7ex = 7br+ = 4/3 and 7bm = 7me = 7dm = 7bm+ = 2/3 (note 
that these cases differ by having two or one points tied to 0) and, as |z| ^ oo 
with I arg z\ < 7r/6, 

(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 

We write the right hand sides as ih^j.{z), . . . ,ih^^_^_{z), and thus these for- 
mulae can be written 

hs{z) ~ ih^iz), (5.9) 

where h'^j.{z) = 2z~^l'^ , ^eJ\z) = 8z, and so on. 

Consider, for simplicity, first the cases H G {br, ex, me, br+} where = 
1/2. We then rewrite <\3.6\\ as, using /5 for the density of \/2 5s. 



hhr{z) 


~ 2iz-V2 




hcx{z) 


~ 8iz 




hhm{z) 


~ 2iTT^/^Z 


-1/4 


^mc(^) 


~ 2i7r^/2^ 


1/4 


hdmiz) 


~ Aiir^/^z 


-1/4 


hhT+{z) 






hhm+{z) 


~ in'/^z- 


1/4 _ 



/7T / Z f'OO 
/ Fo(r,0)e'^«('''^'^'«)drd0 (5.10) 
--it/2 Jo 

where, with 7 = 7=, 

37r-3/2 



Fo{r,e) :-- 



p2ie/3 



^1 

Mr, 0; X, ■■= sec Oe'^ - e'^(e secO/rf/^ - -ir^l^e'^l^. (5.12) 

Remember that ^ is arbitrary; we choose ^ = px^l^ for a positive constant p 
that will be chosen later. Further, make the change of variables r = x^/^s^/'^. 



(sece)V"2/iH(re'^/^), 



(5.11) 
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Thus, 

r'7r/2 noo 



/"TT/Z nOO 

/|(x)=pV/M / Fi(s,6';x)e^'^i(^'^)dsd6' (5.13) 

Je=--K/2 Js=o 



where 



e- x) := —^e^'^/^{sec efs-^/^h^{x^/^s^/^e'^/^), (5.14) 
ipi{s,e) := p(l + i tan 0) -p^/Sg-igie (sec 0)3/2 _^^gie/2^ ^^ -^^^^ 
In particular, 

Reipi{s,9) = p- p^/^s-^ {cos 9)-^/^ --fscos{e/2). (5.16) 
In the cases H G {bm, dm, bm+} when = 1, we obtain similarly 

r7v/2 poo 

/i(x) = p3/22,7/3 / / Fi(s,0;x)e^'^i(^''')dsd0 (5.17) 

Je=-n/2 Js=0 

where 

i7^(s, ^; ^) := J^^e'^/3(sec 9f/h-'/^h^{x'/h^/h''/^) (5.18) 

and is the same as above. 
Consider first 9 = 0; then 

(/7i(s,0) = Rev9i(s,0) = /) - p^/^s"^ - 7S, (5.19) 

which has a maximum at s = sq := p^^^j'^^"^. In order for (so,0) to be a 
saddle point of ipi , we need also 

= ^(.0, 0) = ip - ip'/h-' - lijs = i{p- Ip'/V^') (5.20) 
and thus 

.3y/2.4_ .9^x2_ J9, He{br,ex,br+}, 
^ V2/ V4/ ^9/4, H G {bm,me,dm,bm+}. ^' ^ 

With this choice of p, we find from (j5.15p and (j5.20p that the value at the 
saddle point is 

/ n 3/4 1/2 P [-3, H e {br,ex,br+}, 

o 1—3/4, ^ G |bm, me, dm, bm+|. 

(5.22) 

This yields the constant coefficient in the exponent of the asymptotics. We 
denote this value by —b = —6=, and have thus, using (I5.20p . 

p = 36, p^^/^i/^ ^ 2b. (5.23) 

Further, by (fOHjl and ((5:22]) 

= ,^V.,-./. = 26.- = ll''- ! I '^^'=■'7'; (5.24) 

9/4, G |bm, me, dm, bm+|. 
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The significant part of the integrals in (15.130 and (15.171) comes from the 
square 

Q := {{s,e) : \s - Sol < logx/x, \e\ < logx/x} (5.25) 

around the saddle point, as we will see in Lemma l5.ll below. We consider 
first this square. 

By (I5l4l) . (15181) and ([S^J-dSSl), uniformly for {s,0) G Q, as x ^ oo, 

, , ('i0s^'/\U-^/hl/'), HG{br,ex,me,br+}, 
Fi{s,e;x) = { (5.26) 

[ i+£(l)g-4/3^0 (^4/3^2/3^^ ^ ^ ^^^^ ^^^y 

For the exponential part, we let s = so(l + u/x) and 6 = 2v/x, and note 
that Q corresponds to 

Q' := {{u,v) : |^z| < (logx)/so, 1^1 < (logx)/2}. (5.27) 

A Taylor expansion yields, for {u,v) G Q' , after straightforward computa- 
tions, 

ipi{s,e) = -b-2bu'^x-'^ + 2ibuvx-'^ -bv'^x-'^ + 0{{\u\^ + \v\^)x-^). (5.28) 
Hence, 

^x'Ms,0) ^^^Q ^ 2S0X-2 / / ^-bx^-2bu^+2ibuv-bvHoil) ^^^^ 
Q J Jq' 

j j e-26«^+2i6™-6.2+o(l) ^ ^(1) j 



2sox-^e-''^' ( TT 



25 -16 
-lb b 



-1/2 

+ o{l) 



^2sovr^_2^_,.._ (5.29) 

Further, //Q|e^'^i("'^)| dsd6' is of the same order. Consequently, if we write 

Gi{s,e-x) :=Fi(s,e;x)e^''^i("'^\ 
then ([09]) and (fOGD yield 

[i7|=l5o'^'^/i°(^'/^«o^^)^~'e-^"', {br,ex,me,br+}. 

G,{s,e;x)dsde=i J _ : ^ 

"'"'Q 1 2v/3 7rfc ^0 V )^ ^ ' ^ G {bm, dm,bm-|-}. 

(5.30) 

For the complement := (0, cxd) x (— 7r/2, 7r/2)\(5, we have the following. 
Lemma 5.1. For every N < oo, for large x, 

\Gi{s,d;x)\dsd9 = o(x-^e-^''' 
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We postpone the proof and find from ()5.13p . ()5.17p and ()5.30p . using 

f^s"^/^/ia(x4/3s^/^)x2/3e-^^', E e {br,ex,me,br+}, 
f*{x) ^ J 2v/^ " " ' 

[^s~^^\l{x*/^sl^^)x^/^e-^'=\ E e {bm,dm,bm+}. 

Substituting the functions /i| implicit in (j5.2p - ()5.8p and the values of p, b 
and So given in (j5.2ip - (j5.24p . we finally find 









' 


(5.31) 


/exi 






36\/3 9 
= — ^x^e , 


(5.32) 


f bm ' 




\/3p 1/2 ba;2 




(5.33) 


f* 1 
J me ' 




-1/2 -fcj^a 

~ V o psg 2;e 


3\/3 _3a.2/4 

= -^xe /^ 


(5.34) 


/dm' 






_ 2V3^_32.2/4 


(5.35) 


/br+' 






' 


(5.36) 


/bm+ ' 




\/3/J 1/2 




(5.37) 


^20F^° ^ 


20F 



Recall that these are the densities of V^-Sh. The density of i^s is /h(2^) = 
fi{V^x), and we obtain the leading term of the asymptotics in Theorems 
I1.1H1.71 The leading terms of the asymptotics for ¥{B^ > x) follow by 
integration by parts, as discussed in Section [2l 

It remains to prove Lemma 15.11 We begin by observing that by (j3.7p . 
([5T]) and ([Qjl - dts]) . 

\h{z)\ = 0{\z\ + \z\-'^), |argz|<7r/6. (5.38) 

Hence (j5.14p and (jS.lSp show that, with some margin, 

\Fi{s,e;x)\ <Ci{x'^s~^ + x-'^s-^){cose)-^. (5.39) 

and thus by ()5.16p . for x > 1, 

\Gi{s,e;x)\ < C2x2(cos0)-3(s-l + s-3)e''-'-'^W«"'-'^W^ (5.40) 

where A{9) = p^/'^{cos6)~^/'^ and B{9) := 7cos(6'/2). We integrate over s, 
using the following lemma. 
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Lemma 5.2. Let M > 0. 

(i) If A and B are positive numbers and AB > 1, then 

/ s-''-^e-^'~"-'''ds<C3{M){B/A)^'/\"''^. (5.41) 
Jo 

(ii) If further < 5 < 1, then 

-M-l^-As-^-Bs < C4(M)(i?/yl)*^/2e-(2+^'/2)v^. 



s e 

s-^a/b\ >s^/aJb, s>0 

(5.42) 



Proof, (i): The change of variables s = \J AjB t followed by t i— > t ^ for 
t > 1 yields 

Jo t 

(5.43) 

For t G (g, 1) we write t = 1 — u and use (1 — m)~^ + 1 — u > 2 + u^; hence 
the integral over (|, 1) is bounded by 

C5{M) / e-^(2+«2) ^ C6(M)e-2^^^ 
Jo 

For t G (0, |) we use 

^-M^lg-VABt-V2 < C7(M)(^5)-(*^+l)/2 < C7(M); 

hence the integral over (0, g) is bounded by 

CriM) e-^* /2^t<C7(M)e-3^. (5.44) 
Jo 

(ii): Arguing as in (j5.43p . we see that the integral is bounded by 

/■l/(l+5) , , Af 

The integral over (0, 1/6) is bounded by (j5.44p . and the integral over (1/6, 
1/(1 + 5)) by 

C78(M)e-^(i+^+i/(i+^)) < C8(M)e-^(2+5V2). □ 
Proof of Lemma [3T71 Returning to (|5.40p . we have B{6)/A{6) < ^/ p^^^ and 

A{e)B{e) = p2/3^(cos e)-^/'^ cos(^/2). (5.45) 

Noting that p^/s^ = (25)2 ^^^^ and 
cos z cos a I I cos p 
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we see that 

^A{e)B{e) > 2b + C20^. (5.46) 
Hence Lemma 15.21 applies with A = x'^A{9) and B = x'^B{9) when > 
1/(26) and shows, using (j5.40p . that for every 6 with \0\ < tt/2, 

\Gi{s,e;x)\ds < C9x2(cos^)-=^e^^'-^^^'-2^2x2e2 



= Cgx^{cosey'^e-^''^-^^'=^^\ (5.47) 

For \6\ close to 7r/2, we use instead of (j5.46p 

y^A{e)B(9) > C4(cos 9)-^/\ (5.48) 

another consequence of (j5.45p . Hence, (j5.40p and Lemma l5.2l fi) show that 
if e > is small enough, and \6\ > it/2 — e, then 

\Gi{s,e;x)\ds < Ciox2(cos0)-V^'-'=*^''(™'^^)"'^* 

Moreover, (j5.45p implies that if |^| < 1, say, then 

y/A{e)/B{9) = so + 0{9^). 
Hence, if 16*1 < (logx)/x and \s — so\ > (logx)/x, then, for large x, 

- VmJm > ^ > c,^-^^A{e)/B{e), 

2x X 

and Lemma l5.2f ii) implies, using (j5.46p . that if \9\ < (logx)/x, then 
/ \Gi{s,9;x)\ds < Ci2x2e^^'-^''^'-^''(i°s^)' 

J |s— so|>loga;/a;, s>0 

< Ci3e-^^'-'=^('°s^)'. (5.50) 

The lemma follows by using (|5.50p for 16*1 < (logx)/x, ()5.49p for 16*1 > 7r/2-e, 
and (j5.47p for the remaining 6, and integrating over 9. □ 

6. Higher order terms 

The asymptotics for fs{x) obtained above can be refined to full asymp- 
totic expansions by standard methods and straightforward, but tedious, cal- 
culations. With possible future extensions in view, we find it instructive to 
present two versions of this; the first is more straightforward brute force, 
while the second (in the next section) performs a change of variables leading 
to simpler integrals. 

First, the asymptotics ()4.27p and (j4.28p can be refined into well-known 
asymptotic series [1, 10.4.63,10.4.66 (with a typo in early printings)]; we 
write these as 

Bi'(z) = 7r-V2^V4e2.3/2/3^^(^)^ ^^^2) 
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with 

/3o(^) = 1 + 1^-3/2 

^ 48 4608 



= 1 



48' 



.-3/2 



455 
4608' 



+ ... + 0(^-3iV/2)^ |arg(z)| <7r/3-5, 

(6.3) 

+ ••• + 0(^-3^/2), |arg(z)| <7r/3-(5, 

(6.4) 

where the expansions can be continued to any desired power of z^/'^. 
Similarly, (j4.29p can be refined to an asymptotic series 

BI(z) = 7r-i/2z-3/4e2-'''V3^_^(^)^ (6.5) 

with 

= 1 + f^z-^'^' + + • • • + 0(z-3^/2), I < _ ^. 

(6.6) 

this is easily verified by writing (I4l3l) as Bl{z) = BI(1) + f'^m" {t) dt 
followed by repeated integrations by parts, as in corresponding argument for 
AI(z) in [9|, Appendix A]. The coefficients in (16. 6p are easily found noting 
that a formal differentiation of (|6.5p yields (j6.ip . (They are the numbers 
denoted in [o'].) 

Hence, by (j4.30p - (|4.36p . (|5.9p can be refined to, for | argzj < 7r/6, 



h^{z) = \hl{z) (hl{z)+o( e-^-^'"/^) ) , 



(6.7) 



where 



hl^iz) 


:=/3i(^)-' 


= 1 + 


7 

24^" 


-3/2 


+ ..., 


(6.8) 




:=/3i(z)/3o(^)~' 


= 1 - 


11 

24^" 


-3/2 


+ ..., 


(6.9) 


^bm(^) 


:=/3i(^)-' 


= 1 + 


7 

48^" 


-3/2 


+ ..., 


(6.10) 




:= /3o(^)~' 


= 1 - 


5 

— z 
48 


-3/2 


+ 


(6.11) 


^dm(^) 


:= [i-Mz)-^ 


= 1 + 


31 

48^" 


-3/2 


+ ..., 


(6.12) 




:= {mz)+l3^{z))/2) 


-' = 1 + 


1 

24^" 


-3/2 


+ ..., 


(6.13) 


^bm+(^) 


:= {mz)+Pi{z))/2) 


-' = 1 + 


1 

48^" 


-3/2 


+ ... 


(6.14) 



By (16. 3p and (16. 4p . /i4 has an asymptotic series expansion 

hl{z) = 1 + dfz-'^/^ + dfz-^ + ■■■ + 0{z-^^/^), I arg(z)| < tt/3 - 6, 

(6.15) 

for some readily computed coefficients d^; moreover, we can clearly ignore 
the O term in (16.71). 
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Next, by Lemma [5?T| it suffices to consider (s, 9) £ Q in ()5.13p and ()5.17p . 
We use (j6.7p and (j6.15p in ()5.14p and (|5.18p and obtain, for example, for 
S = ex, 

Fi e;x) = ^ (e*^(sec dfs~' - ^e''/' (sec 9fs-'x-' + ■■■ + ©(x'^^)) . 

(6.16) 

We substitute s = sq{1 + u/x) and 9 = 2v/x as above and obtain by Taylor 
expansions a series in x~^ (with a prefactor x^^^) where the coefficients are 
polynomials in u and v. 

Similarly, the Taylor expansion (|5.28p can be continued; we write the 
remainder term as i?(u, v] x) and have 

R{u, v; x) = rsiu, + • • • + ©(x'^^-^), (6.17) 

for some polynomials ri.{u,v). Another Taylor expansion then yields 

^x^R{u,v;x) ^ ^ . . . ^ 0(x-2^), (6.18) 

for some polynomials r'^{u,v). We multiply this, the expansion of Fi and 
the main term exp(— — 2bu^ + 2ibuv — bv'^), and integrate over Q; we 
may extend the integration domain to with a negligible error. This yields 
an asymptotic expansion for f^{x), and thus for f^ix), where the leading 
term found above is multiplied by a series in x~^, up to any desired power. 
Furthermore, it is easily seen that all coefficients for odd powers of x~^ 
vanish, since they are given by the integrals of an odd functions of u and v; 
hence this is really an asymptotic series in 

We obtain the explicit expansions for /h(2;) in Theorems ll.lffTTTl bv cal- 
culations with Maple. The asymptotics for ¥{Bs > x) follow by integration 
by parts, see Section [2j 

Remark 6.1. In particular, since h^j.j^ = ^/ibj. and hY^^_^_ = ^/ibj^) the 
leading terms for br+ and bm+ differ from those of br and bm by a factor ^ 
as discussed in Remark 11.81 The second order terms in are different, as is 
seen above; more precisely, = /ij^j, — jz~^^'^ + 0{z^^) and hl^^^ = /ibm ~ 
12,-3/2 _|_ (9(_2~^); it is easily seen that if we ignore terms beyond the second, 
this difference transfers into factors 1 — (4so)~^x~^ and 1 — (8so)~^x~^, 
respectively, for /5, which in both cases equals 1 — j^x'"^, and thus a factor 
1 — ^x~^ for /s, which explains the difference between the second order 
terms in /br or /bm and 2/br+ or 2/bm+; cf. again Remark 11.81 

7. Higher order terms, version II 

Our second version of the saddle point method leads to simpler calcu- 
lations (see for instance, Bleistein and Handelsman [2]). We illustrate it 
with Bex] the other Brownian areas are treated similarly. We use again 
(I5.13p . and recall that by Lemma 15.1^ it suffices to consider {s,9) close to 
(so,0) = (9/2,0). 
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We make first the substitution s = |(sec^)'^/^u ^ (this is not necessary, 
but makes the integral more similar to Tolmatz' versions). This transforms 
(f5l3]l into 

/'7r/2 /"OO 

fUx)= / F2{u,e-x)e-''P^^-^'Uude, (7.1) 

Je=--K/2 Ju=o 

where, by (KIM . (KIM . dOTD . (K7\\ . ([53]), (US]), for u bounded, at least, 

~ 7r3/2n(cOS0)3 " 2^3/2(cos^)3/2 +^^^^' ^^"^^ 

^2(n,0) = 9(1 + itan^) - 6txe'^ - + itang)3/2 ^ ^^^^^ 

u 

The saddle point is now (u, 9) = (1, 0), and in a neighbourhood we have, cf. 
(IS2SD, with v = u-l, 

ip2{u, 9) = -3- - 3w9 - + 0{\vf + |0|3). (7.5) 

The function ip2 has a non-degenerate critical point at (1,0), and by the 
Morse lemma, see e.g. Milnor |16l . Lemma 2.2], we can make a complex 
analytic change of variables in a neighbourhood of (1,0) such that in the 
new variables (/?2 + 3 becomes a diagonal quadratic form. (The Morse lemma 
is usually stated for real variables, but the standard proof in e.g. [l^ applies 
to the complex case too.) The quadratic part of (17. 5p is diagonalized by 
{v, 9) with V = V — \9 /A] we may thus choose the new variables u and 9 such 
that u ^ V and 9 ^ 9 &t the critical point, and thus 

u = l+u-\9/A + 0{\u\^ + \9\^), (7.6) 
= ^ + O(|n|2 + |^|2)^ (7.7) 

(^2(n,e) = -3-6u2-|^^ (7.8) 

Note that the new coordinates are not uniquely determined; we will later 
use this and simplify by letting some Taylor coefficients be 0. In the new 
coordinates, (|7.ip yields, 

/ex(^) ~ / / F3(n,^~;x)e-3-'-^'(6"'+l^'')j(n,0~)dude; (7.9) 



where F3 is obtained by substituting u = u{u, 9) and 9 = 9{u, 9) in (j7.2p 
and J{u,9) = f||| — is the Jacobian. Recall that, up to a neghgible 

error, we only have to integrate in (j7.ip over a small disc, say with radius 
logx/x; this becomes in the new coordinates a surface in as the integra- 
tion domain in (|7.9p . The next step is to replace this integration domain 
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by, for example, the disc {{u,6) £ : |up + \6\'^ < (logx/x)^}, in anafogy 
with the much more standard change of integration contour in one complex 
variable. To verify the change of integration domain, note that if F(zi,Z2) 
is any analytic function of two complex variables, then F{zi, Z2) dzi A dz2 
is a closed differential form in (regarded as a real manifold of dimension 
four), and thus Jg^^ 22) dzi A dz2 = by Stokes' theorem for any com- 
pact submanifold M with boundary dM. In our case, it follows that the 
difference between the integrals over the two domains equals an integral over 
boundary terms at a distance ^ logx/x from the origin, which is negligible. 
(The careful reader may parametrize the two domains by suitable mappings 
ipo,jpi : U ^ C'^ , where U is the unit disc in M?, and apply Stokes' theorem 
to the cylinder U x [0, 1] and the pullback of F{zi, Z2) dzi A dz2 by the map 
{w,t)^ {l-t)i;Q{w)+t^i{w).) 

We next change variable again io w = xu, t = x9, and obtain by (|7.9p 

ft,{x) ~ x-V^^' j j F3{w/x,t/x;x)J{w/x,t/x)e-'^'"'-'i'' dwdt, (7.10) 

integrating over {w,t) € with, say, w'^ + t'^ < (logx)'^. To obtain the 
desired asymptotics for f*^, and thus for /ex, we mechanically expand F3 
and J in Taylor series up to any desired order and compute the resulting 
Gaussian integrals, extending the integration domains to M^. 

We illustrate this by giving the details for the first two terms in (II. ip . We 
have, cf. (j7.6p and (|7.7p . expansions 

u = l + u- i(9/4 + aiu^ + a2u0 + 0{\uf + |^|^), 

e = 9 + 03^^ + a^uO + 0{\uf + \9f), 

where we, as we may, have chosen two Taylor coefficients to be 0. To de- 
termine ai, . . . , 04, we substitute into f2iu, 0). We obtain from ()7.4p . up to 
terms of order three, 

ip2{u, 9) 3 - [6^2 + 961^8] + (6 - I2ai)u^ + (-3ia4 - 12a2 - lb\/2)9u^ 

( 9q4 . 33\ ~2 / 9a3 151 \ ~n 

Annihilating the coefficients, cf. (j7.8p . leads to a linear system, the solution 
of which is 

ai = 1/2, 02 = -831/72, 03 = 51/24, = 19/9. 

This leads to the Jacobian 

Furthermore, by (j7.3p . with w = xu and t = x9, 

^ / ;t N ^/ ^ N 162x'^ 162x^/5. \ ^/ 9, 9 9x\ 
F3iu,9;x) = F2iu,9;x) + —^(^-it - wj + 0{x\l + + t^)) . 
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Integrating in ()7.10p yields the leading term 

36^3 2„-3x2 



/ex(^)~7lk^^ (7.11) 



together with correction terms of order xe~^^ that all vanish by symmetry, 
since they involve integrals of odd functions, plus a remainder term of order 

g-3a; _ 

n 2 

The next term in the expansion of e f*^ is thus the constant term. To 
find it, we try, again setting some Taylor coefficients to as we may, 

n ~ 1 + (il - ie/A) + u{aiu + a2d) + u{Piu^ + I32u9 + PsO"^), 
We obtain now 

ip2{u,e) 3- [6^^96*78] + (3/2- 12/3i)u'^ + (-131i/6-3i/34- 12/32)^^61 

+ (-9/34/4 + 2627/288 - 12/33 - 3il35)u'^e'^ + (-9/^5/4 + 5351/96 - 3i(36)u9^ 

+ (-1283/768 - 9/36/4)^^. 

We set for instance /34 = 0. This gives 

/3i = 1/8, /32 = -131i/72,/33 = 16867/10368, /^s = 44931/1296, /Je = -1283/1728. 

The Jacobian becomes 

. ~, f 28u 5ie\ /179 9 2405. ~ 379 

J(u, 6')~1+ + + iuO 

^ ' ^ V ^ 24 y V 72 648 384 

The first term in (j7.3p becomes 

162j;'^ 162x^/5.- .\ / 2 1143. . 621 -3 



7r3/2 I' 



+ -^HTT {-:iO-u] + — - 81«^ + -—iu9 + 



7r3/2 7r3/2 V4 / 7r3/2V 4 
and the second is 

33x2 



2^3/2 • 

Collecting terms, the coefficient of x^ in F-i{w/x^t/x)J{w/x,t/x) equals 
1296^2 - 99248iwt + 2565*^ 33 



647r3/2 27r3/2 

8 

„2 



Multiplying by exp(— 611;^ ~ |*^) ^'^'^ integrating yields the contribution 



to the integral in (j7.10p . So, finally, combining (j7.1ip and ()7.12p . 



3l/2g-3x2 

fcAx) — [36x 



which fits with the first two terms for fexi^) in Theorem II. li More terms 
can be found in a mechanical way. 
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8. Proof of Theorem 13.11 

Let T ~ r(z/) be a Gamma distributed random variable independent of 
X and let Xt := T^/^X. Then T has the density Tiiy^H^-^e-^ t > 0, and 
thus Xt has, using (13. the Laplace transform 

/■OO 

Jo 

= ^-^3^(^-2/3^^ n>0. (8.1) 



By ()8.ip and our assumption on ipT extends to an analytic function in 
C \ (oo,0]. Furthermore, Xt has a density g on (0,oo), because r3/2 has, 
and it is easily verified that this density is continuous. We next use Laplace 
inversion for Xt- The Laplace transform ipT is, by ()8.ip . not absolutely 
integrable on vertical lines (at least not in our cases, where ^'(2:) is bounded 
away from as z — > 0), so we will use the following form of the Laplace 
inversion formula, assuming only conditional convergence of the integral. 

Lemma 8.1. Let h be a measurable function on M. Suppose that the Laplace 
transform h{z) := h{y)e''^y dy exists in a strip a < Re z < b, and that 

a G (a, b) is a real number such that the generalized integral f^^^^ e^^h{z) dz 
exists in the sense that the limit limy^^oo sa exists, where 



t-cr+iA 

SA ■= / e''''h{z) dz. 

Ja-iA 



If further x is a continuity point (or, more generally, a Lebesgue point) of 
h, then 

/ e''''h(z) dz := lim = 27ri/i(x). 

Ja-ioo 

Proof. By considering instead e~'^^/i(y), we may suppose that a = 0. In this 
case, h is integrable and h{\t) = h(t), the Fourier transform of h, and the 
result is a classical result on Fourier inversion. (It is the analogue for Fourier 
transforms of the more well-known fact that if a Fourier series converges at 
a continuity (or Lebesgue) point of the function, then the limit equals the 
function value.) For a proof, note that if sa converges as ^ ^ 00, then so 
does the Abel mean ye~'^^SA dA as y ^ 0, and this Abel mean equals 27ri 
times the Poisson integral iT~^y{u'^ + y'^)~~^h{x — u) du, which converges 
to h{x). □ 
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We verify the condition of the lemma with h = g and a = 1, recalling 
that g = ^T- Thus, by I^A\i . 

/•l+iA rl+iA 
Jl-iA Jl-iA 

We may here change the integration path from the straight line segment 
[1 — iA, 1 + iA] to the path consisting of the following seven parts: 

7i: the line segment [1 — iA, — A — iA], 

721 the line segment [—A — iA, — A — iO], 

73: the line segment [—A — iO, — e — iO], 
74: the circle {ee'* : t £ [— 7r,7r]}. 

75: the line segment [—e + iO, — A + iO], 

76: the line segment [—A + iO, — A + iA], 

77: the line segment [—A + iA, 1 + iA]. 

(Here, 73 could formally be interpreted as the line segment [—A—irj, — y^e^ — r/^- 
irj] for a small positive t], taking the limit of the integral as r/ — > 0, and 
similarly for the other parts with ibiO.) Letting A 00, we see that we 
essentially change the integration path from a vertical line to a Hankel con- 
tour; however, we do this carefully since, as said above, the integral along 
the vertical line is not absolutely convergent. 
We now first let e ^ 0. By ([OH . 



,x.^-2v/3^(^-2/3) ^^^Q ^ 0, 

'74 

and, again by ()3.3p . the integrals along 73 and 75 converge to the absolutely 
convergent integrals 

J-A-iO Jo 

and 

JiO Jo 

which together make 

l-A 
Jo 

Hence, for every A > 0, 

sa = Ia+([ +[ +[ + [)e^^z'"'/'^{z~'/')dz. 
Now let A ^ 00. By i^M . 

j gX^^-2,./3^(^-2/3) d^ ^ ^ (^j^ ^xt ^ ^(^^^ 
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and similarly /^^ = o(l), /^^ = o(l), f^^ = o(l). Finally, Ia loo, and 
Lemma |8. II applies and yields the following. 

Lemma 8.2. For every x > 0, we have 

1 

9{X) = ^ g-xp^-2./3^* (^-2/3) (8.2) 

where the integral is absolutely convergent by ()3.7p and (jS.Sp . 
By the change of variables p = u~^^'^, (j8.2p may be rewritten as 

q /"OO 

g(x) = — e-^""'^'n'^-5/2^*(u)du, x>0. (8.3) 
4vri Jo 

We can here, using ()3.7p and (j3.8p . change the integration path from the 
positive real axis to the line {re"'' : r > 0}, for every fixed with \ip\ < ^. 
Consequently, we further have, for x > and \ip\ < ^, 

q fCO 

g(x) = —e^''-^'^y'^ exp(-e-3'^/2^r-3/2)r''-5/2^*(re''^)dr. (8.4) 
47ri 7o 

The right hand side of (|8.4p is an analytic function of x in the sector {x : 
I argx — 3(/?/2| < vr/2}, which contains the positive real axis; together, these 
thus define an analytic extension of g{x) to the sector | argx| < 37r/4 such 
that (j8.4p holds whenever | argx| < 37r/4, \ip\ < ^ and | argx — 3(/?/2| < tt/2. 

We next find the density of X from g by another Laplace inversion. As- 
sume first, for simplicity, that we already know that X has a continuous 
density / on (0,oo). Then t^^'^X has the density t~'^^'^ /{t^'^^'^x), and thus 
(using t = x^/^s), for x > 0, 

/•oo 

g{x) = Tiu)-' / f'-'e-H-^/^fit-^/^x) dt 
Jo 

/•oo 

= Tii^r'x^''/'^-' / e-^'''^s'^-5/2/(,-3/2) d,. (8.5) 
Jo 

Let 

F{s) := s^-^/^f{s-^/^). (8.6) 
Then (18. 5|) can be written, with x = y'^^'^, 

/•oo 

g{y^'^) = T[v)-^y'^-'^'^ / e-y'F{s) ds, y > 0. (8.7) 

In other words, F has the Laplace transform 

F(y) := / e-y'F{s) ds = T{v)y^l^-^ g{y^'^), y > 0. (8.8) 
JO 

Since this is finite for all y > 0, the Laplace transform F is analytic in the 
half-plane Rey > 0. Hence, using our analytic exytension oi g io \ argz| < 
37r/4,, (18. Sp holds for all y with Rey > 0. Consequently, by standard Laplace 
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inversion, for every s > and every > such that the integrals are 
absolutely (or even conditionally, see Lemma IS.ip convergent, 

F{s) = — / e^ypiy) dy = / e^^y'/'^^ g{y^/^) dy. (8.9) 

We have for (mainly notational) simplicity assumed that X has a density. In 
general, we may replace the density function / in ()8.5|) and ()8.6p by a prob- 
ability measure /i (with suitable interpretations; we identify here absolutely 
continuous measures and their densities as in the theory of distributions). 
Then F is a (positive) measure on (0, oo), and its Laplace transform is still 
given by (j8.8p . The fact, proved below, that F is absolutely integrable on a 
vertical line Rey = ^ implies by standard Fourier analysis that F actually 
is the continuous function given by (j8.9p . and thus the measure n too is a 
continuous function; i.e., X has a continuous density / as asserted, and (|8.5p 
and hold. 

We change variables in (jS.Op to 9 := argy G (— 7r/2, 7r/2), that is y = 
^(l + itan6') = ^sec(6')e'^ We further express ff(y^/^) by ([SID with (/? = 61/3 
(which satisfies the conditions above for ()8.4p l: this yields, assuming absolute 
convergence of the double integral, 

F{s) = ^ r expfCs(l + itan^) - e'' {Csec9f/\-'/^) 

g(l-2./3)ie^5/2-.(gg^^)7/2-.^.-5/2^*(^gie/3) (g^^^p) 

To verify absolute convergence of this double integral, take absolute values 
inside the integral. Since 

Re(e'^(^sec0)3/V"3/2) > fl\^^^e)'l\-^l\ 

the resulting integral is, using (j3.5p and (j3.7p . for fixed 5 and ^ bounded by 

C7i4(s,0 / / e-«'^'(^^^^)'^'^"'^'(sece)^/2-.^.-5/2^.^(^-.^^-6)^^^^^ 
Jo JO 

We split this double integral into the two parts: < < 1 and 1 < < oo. 
For < < 1, sec 6* is bounded above and below, and it is easy to see that 
the integral is finite. For > 1, tan^ < sec0 < 2tan0, and with t = tan9 
we obtain at most 

/oo poo 

Substituting t = r^u, we find that this is at most 

/"CXD /'OO 

Ci5 / e-«'^'"'^'n3/2-- dTX / r5-2- min(r-^ r^^) dr < oo. 
Jo Jo 

This verifies absolute convergence of the double integral in (]8.10p for every 
^ > 0, which implies absolute convergence of the integrals in ()8.9p . Conse- 
quently, (18. 9p and (j8.10p are valid for every s > and ^ > 0. We now put 
s = x"^/^ in (fSTTOl) and obtain by ([8^ the sought result ([3^ . 
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Remark 8.3. We have chosen ip = 9/3 , which leads to (13. 6p and, see Re- 



mark 13.31 the formulas by Tolmatz [23, l2J, |25| . Other choices of (p are 



possible and lead to variations of the inversion formula ()3.6[) . In particular, 
it may be noted that we may take p = for, say, \6\ < 1; this yields a 
formula that, apart from a small contribution for \6\ > 7r/4, involves ^*{x) 
for real x only. However, we do not find that this or any other variation 
of (|3.6p simplifies the application of the saddle method, and we leave these 
versions to the interested reader. 

9. Moment asymptotics 

Suppose that X is a positive random variable with a density function / 
satisfying (|2.ip . Then, as r ^ oo, using Stirling's formula, 

oo 

EX"~ / ax'+^e"''^' dx 







oo 



y{r+a+l)/2-l^-by 







= "^-(r+a+l)/2p/ ^ + Q + l 

2 V 2 

= aV^(26)-(-+i)/2r°/2f^Y'^'. (9.1) 

\2eb/ 

(It is easily seen, by an integration by parts, that the same result follows 
from the weaker assumption ()2.2p .) 

For the Brownian areas studied in this paper. Theorems I1.1H1.7I thus 
imply the following. 

Corollary 9.1. As n ^ oo, 

EB2^^3V2n[—) , (9.2) 

EBli,^V2[—) , (9.3) 

EBl^,,r.V2(^Y^\ (9.4) 



3e 



EB^^r^V3^n'/'[^f\ (9.5) 
EeV~2^ (^Y^\ (9.6) 



V3e. 

EB^^^^^i-^) ' , (9.7) 



1 / n \ "/2 
72^12^^ 
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Most of these results have been found ear her : (|9.2I) by Takacs UM, (19-30 by 
Takacs and Tolmatz (2i| JO]) by Takacs [2l|^d Tolmatz 0, ([93]) by 



Takacs [22] , ()9.6p by Janson [9] , (19. 7[) by Tolmatz [251] ; Takacs used recursion 
formulas derived by other methods, while Tolmatz used the method followed 
here. Note that, as r emarked by Tolmatz [2I, IESbr+ ~ ^ ^^-^br and similarly 
^>3l^m+ ~ llE^^^, cf. Remark [13 

In the opposite direction, we do not know any way to get precise asymp- 
totics of the form (|2.ip or (j2.2p from moment asymptotics, but, as observed 
by Csorgo, Shi and Yor the much weaker estimate (|1.3p and its analogue 
for other Brownian areas can be obtained by the following special case of 
results by Davies [3] and Kasahara [HJ. (See [sl. Theorem 4.5] for a more 
general version with an arbitrary power instead of in the exponent.) 

Proposition 9.2. If X is a positive random variable and b > 0, then the 
following are equivalent: 

-lnP(X > x) ~ X ^ 00, 



ln(Ee*^)~4^^ t ^ 00. 

Returning to (j9.ip . we obtain in the same way more precise asymptotics 
for the moments if we are given an asymptotic series for / or ¥{X > x). 
For simplicity, we consider only the next term, but the calculations can be 
extended to an asymptotic expansion with any number of terms. Thus, 
suppose that, as for the Brownian areas, (j2.ip is sharpened to (j2.3p with 
N > 2. Then, also using further terms in Stirling's formula, 



0(5-n/2r/- + «-3 



l6-(n+a-M)/2r(!i + ^) ( + ,^5! + Q 

2 V 2 / \ n 



in-') 



27r(26)"("+i)/2n"/2(^)"^' 



[oq + ^ao " ^ ^ + y + n ^ + 0{n 



2^ 



In particular, for the Brownian excursion, where by Theorem ll.il (j2.3p holds 
with a = 2, 6 = 6, Oo = l^^J^j'K and 02 = — S-y/G/vr, 
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If we, following Takacs [19(], introduce Kn defined by 

_ 4V?2-/^n! 
r((3n-l)/2) 

further applications of Stirling's formula shows that (j9.9p is equivalent to 
A-„M2.)-"VV^(|)"(l-^.0(O), ,8.10) 



Again, the leading term is given by Takacs [19|, in the equivalent form 

1 /3\" 

7 {n-iy. asn^cx). (9.11) 



27r V4 

Takacs (l9| further gave the recursion formula (with Kq = —1/2) 

3n - 4 ""^ 
Kn = + KjKn-j, n > 1, (9.12) 

i=i 

It is easy to obtain from (j9.1ip and (j9.12p the refined asymptotics 

^" = 2^(i)"("-i>'('-is;+°<""')' 

which is equivalent to (jD.lOp . and, by recursion, ()9.13p can be extended 
to an asymptotic expansion of arbitrary length. (Another method to ob- 
tain an asymptotic expansion of is given by Kearney, Majumdar and 
Martin [12] •) Hence (j9.9p (also with further terms) can, alternatively, be 
derived from (j9.1ip and (j9.12p by straightforward calculations. However, 
as said above, we do not know any way to derive Theorem 11.11 from this. 
(Nevertheless, the calculations above serve as a check of the coefficients in 
Theoremll.il) 
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Appendix A. Proof that ^/zk\{z) - A.]'{z) has no zeros 

The double Laplace transforms for the positive part areas ;Bbr+ and Sbm+ 
have both the denominator y^Ai(2;) — Ai'(z), and it is important that this 
function has no zeros, whence ^'br+ and ^j!!, _|_ are analytic in the slit plane 
C \ (—00,0]. This was proved by Tolmatz [25| (for the same reason), but 
we give here an alternative proof that does not need the careful numerical 
integration done by Tolmatz. (Our proof is, like Tolmatz', based on the 
argument principle.) 

Lemma A.l (Tolmatz (25I]). The function y^Ai(2;) — Ai'(2;) is non-zero for 
all z = re'^ with r > and \9\ < vr. 
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Proof. We use the notations 

fiz) 
9{z) 



2 3/2 
3^ , 



vr (ViAi(z) - Ai'(z)) , 



Note that these functions are analytic in the sht plane C \ (— oo,0] and 
extend continuously to (— oo,0] from each side, so we can regard them as 
continuous functions of re'^ with r > and — vr < r < vr, where we regard the 
two sides re^"^ = — r it iO of the negative real axis as different. (The reader 
that dislikes this can reformulate the proof and study zAi(z^) — Ai'(z^) for 
Rez > 0; this avoids the ambiguities of square roots.) 

We will use the argument principle on g{z) and the contour consisting 
of the interval from to — — 10 along the lower side of the negative real 
axis, the circle i?e'^ for — vr < < ir and the interval from —R + iO back to 
0, where i? is a large real number. 

First, fix a small 5 > 0. By (jiTMl) and (f05|) . as \z\ 00, 



arg z\ < TT — 6. 



Next, assume < argz < 27r/3— 5. Note that then arg(— 2) 
{—TT, — 7r/3 — 5) and thus 



a/2 



-IZ 



'i\ c{-z) = e-('/'y-az) = iaz). 



(A.l) 

arg(z) — TT G 

(A.2) 



00 with I arg 2; I < 2it/3 — 6 the expansions 



Furthermore, we have as \z\ 
Q, 10.4.60,10.4.62] 

Ai(-z) = ^-1/2,-1/4 (^sm{C{z) + ^) (1 + OiC')) - cos(C(^) + ^) • 0{C') 

cosiCiz) + 7) (1 + OiC')) + cos(C(z) + j) • OiC' 



-1/2^1/4 



Ai'( — z) = TT z 

and thus 



fi-z) 



.1/4 



cos(C(z) + 7) (1 + 0{C')) - isin(C(z) + 7) (1 + 0{C' 



In the range arg 2; G (0,27r/3 — 5), further 9C(-2) > 0, and thus by Euler's 
formulas 



cos(C(z) + 7) + sin(C(z) + 7) 



< 



< 2e^«^) 



-z)V4e-C(-^). 



Hence, as \z\ — > 00 with argz G (0, 27r/3 — 5), using ()A.2p . 
f(-z) = ^i/4e-i(CW+-/4) (1 + oiC')) ^ 

Consequently, (jA.lj) holds as \z\ — > cxd with — vr < argz < — 7r/3 — 6 too. 
Since /(z) = f{z), it holds for vr/S + 5 < argz < vr too, and combining the 
three ranges, we see that as \z\ 00, for all | argz| < vr, 

/(z)~zi/V^(^), (A.3) 

and thus 



9{z) 



.1/4, 



1 + 0(1)). 



(A.4) 
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Consider now f(z) on the lower side of the negative real axis, i.e. for 
z = re~"^ = — r — iO, r > 0. Note that then Ai(z) and Ai'(z) are real 
and z^/^ purely imaginary. Since Ai and Ai' have no common zeros, and 
Ai'(O) 7^ 0, /(— r — iO) 7^ 0. Moreover, /(O) > 0, and as r grows from to 00, 
/(— r — iO) is real at r = and at the zeros — r = of Ai, and imaginary at 
the zeros — r = of Ai'. Consider continuous determinations of arg f{z) and 
argg{z) along the neagtive real axis, starting with arg/(0) = argg{z) = 0. 
It is easily seen that arg/(2;) then is — 7r/2 for z = a'^, —tt for z = ai, and 
so on, with arg/(a/c — 10) = —kir. Furthermore, for z = —r — iO, 

argg(z) = arg/(z) + 9C(^) = arg/(z) + §9^3/2 ^ ^^g^(^) ^ 2^3/2_ 
In particular, using the asymptotic formula [l|, 10.4.94] for the Airy zeros 

argg(afc) = -T^k + fjafcl^/^ = -vr/j + ^ 37r(4fe - 1) ^ 0{k~^)) 

= _^ + 0(fe-i). (A.5) 

Consider now a continuous determination of &Tgg{z) along the contour 
with R = \ak\ for a large k. On the part from to —i? — 10, the argument 
decreases by — 7r/4 + 0{k~^) by (jA.Sp . and on the half-circle from —R — iO 
to R, it increases by (|A.4p by 7r/4 + o(l), so the total change from to i? is 
0(1), i.e., tends to as A: —> 00. Since furthermore g{R) > 0, the change is a 
multiple of 27r, and thus exactly for large k. By symmetry, the change of 
the argument on the remaining half of 7^ is the same, so the total change 
along 7/j is 0, which proves that g{z) has no zero inside 7/j for R = |afc| with 
k large. Letting k — > 00, we see that g{z) has no zeros. □ 
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